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Abstract. Various types of metal-insulator transitions are discussed to find condi- 
tions for which an ideal surface of a bulk insulator is metallic. It is argued that for 
the correlation- driven Mott metal-insulator transition the surface phase diagram 
should be expected to have the same topology as the phase diagram for magnetic 
order at surfaces: The corresponding linearized mean-field descriptions, a simplified 
dynamical mean-field theory of the Hubbard model and the Weiss mean-field the- 
ory for the Ising model, are found to be formally equivalent. A new kind of surface 
state appears in the low-energy part of the one-particle excitation spectrum as a 
precursor effect of the Mott transition. 



The Mott metal-insulator transition at a crystal surface is a subject that 
touches different areas in solid-state theory which arc usually treated as being 
disjoined: many-body theory of the correlation-driven metal-insulator tran- 
sition, the general theory of surface phase transitions, and the theory of 
electronic surface states. It is the intention of the present paper to show that 
a corresponding combination of different concepts can be fruitful and allows 
for some new theoretical predictions. 

1 Surface Phase Transitions 

The large variety of novel and interesting phenomena in surface physics is 
closely related with the occurrence of surface phase transitions. As has been 
pointed out by Mills , the surface of a system may undergo a phase tran- 
sition at a critical temperature T CjS being different from the bulk critical 
temperature T c , i.e. the surface may undergo its own phase transition. Criti- 
cal exponents, for example, can be defined and determined which are specific 
for the transition at the surface and which cannot be fully reduced to the 
bulk critical exponents (|]||. Different kinds of surface phase transitions are 
conceivable and have been found, e.g. structural transitions, such as deviating 
geometrical order of the atoms near the surface of a single crystal (surface 
reconstruction) , the loss of long-range crystalline order at the surface prior to 
a bulk melting transition (surface melting) or enrichment of one component 
at the surface of a solid binary alloy (surface segregation) M . Typical exam- 
ples for surface phase transitions are also found among magnetically ordered 
systems: For example, the (0001) surface of ferromagnetic Gd is believed to 
have a Curie temperature which is higher than the bulk Tq p|. 
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Different types of surface phase transitions can be described in a qual- 
itative but consistent way by means of classical Landau theory ]||| or by 
lattice mean-field approaches which may be considered as coarse-grained real- 
izations of the Landau theory. Especially, mean-field approaches to localized- 
spin models, such as the Ising or Heisenberg model, are frequently consid- 
ered in this context j^]. For surface geometries there are a number of non- 
trivial results predicted by Landau or mean-field theory, such as temperature- 
dependent order-parameter profiles, which may give a surprisingly good de- 
scription of experimental data (see Ref. @, for example). 

Typically, the surface undergoes the phase transition at the same temper- 
ature as in the bulk, T CjS = T c , if the local (structural, electronic, magnetic) 
environment remains unchanged, while T CjS > T c if there is a perturbation A 
at the surface exceeding a certain critical value A c . For A > A c and temper- 
atures T c < T < T CjS there is an ordered D — 1-dimensional surface coexisting 
with a disordered /^-dimensional bulk. More complicated phase diagrams are 
obtained in the case of multi-critical behavior, e.g. when the long-range or- 
der at the surface has a character different from the long-range order in the 
bulk (surface reconstruction, anti-ferromagnetic surface of a ferromagnetic 
bulk, etc). The Landau T-A phase diagram should be qualitatively correct 
whenever the D — 1-dimensional system can support independent order |^,[^] . 

2 Metal-Insulator Transitions 

The concept of a surface phase transition and the corresponding Landau 
theory seems to extend straightforwardly to a certain kind of metal-insulator 
transitions, namely those which accompany an order-disorder thermodynamic 
phase transition (see Ref. || for an overview). The thermodynamic phase 
transition will be considered at the T = quantum-critical point as due 
to ubiquitous thermal activation processes, a strict definition of a metal- 
insulator transition is possible at zero temperature only. It is well known that 
the formation of an ordered state may result in a gap for charge excitations 
as, for example, in the Peierls transition or in the transition to an anti- 
ferromagnetic state: Consider the typical example of a bipartite lattice and 
two-sublattice long-range order causing a doubling of the unit cell. For a non- 
degenerate band the splitting at the boundary of the new Brillouin zone will 
lead to a gap and, in the case of half-filling, to an insulating ground state. 

Now, for sublattice order at the surface of a disordered bulk one would 
have the (naive) expectation that an insulating surface could coexist with 
a metallic bulk. This, however, is clearly impossible as a finite bulk density 
of states at the Fermi energy will always induce a non-zero, though possibly 
low density of states in the surface region. Likewise, the reverse scenario 
is impossible either: Namely, to realize a metallic surface phase of a bulk 
insulator caused by a thermodynamic phase transition, a disordered surface 
would have to coexist with an ordered bulk which, in general, is forbidden by 
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strict arguments (though a magnetic "dead-layer" scenario is found under 
somewhat exotic circumstances in a D = 2 q-state surface Potts model 0). 

Besides the thermodynamic transitions, there is a second important class 
of metal-insulator transitions, namely quantum-phase transitions H. Essen- 
tially these take place at T = only and are not associated with any symme- 
try breaking. Important examples are the transition from a metal to a nor- 
mal band insulator and the Mott-Hubbard transition from a metallic Fermi 
liquid to a Mott insulator. While the former can be understood within an 
independent-electron model, correlation effects are constitutive for the lat- 
ter. Referring to a quantum-phase metal-insulator transition, it is very well 
feasible that the surface of a bulk insulator is metallic. 

3 Surface States 

Fig. [I] shows a possible electron density of states for this situation. To have 
a metallic surface of an insulator, the density of states at the Fermi energy 
E-p must be finite at the surface while, with respect to the bulk states, Ep 
should lie within a band gap. Note that this necessarily implies the existence 
of a partially filled surface state: The appearance of a surface state at the 
Fermi energy is crucial to get a metallic surface of an insulator. 

Two possible origins of electronic surface states are well known [jTo| : (i) 
Image-potential states may arise as Rydberg-like states in the long-ranged 
— I j Az image potential which is due to the polarization charge that is induced 
by an electron approaching a conductive surface. The electron can be trapped 
between the image-potential surface barrier and the bulk barrier due to a bulk 
band gap. (ii) Crystal-induced surface states originate from the mere crystal 
termination. For an ideal unreconstructed crystal surface a further distinction 
between Tamm and Shocklcy states is meaningful: Shockley states appear 
within a hybridizational band gap which may open when the boundaries of 
two bulk bands have crossed as a function of decreasing lattice constant. 
Tamm-like surface states are due to the surface change of the one-electron 
potential and always lie near the bulk band from which they originate. 

In fact, a crystal-induced surface state may give rise to a metallic surface 
phase as is demonstrated by the following two examples: Due the reconstruc- 
tion of the Si(lll)-(7 x 7) surface there is a partially occupied surface state 
which is consistent with the observed metallic behavior for this surface (cf. 
the discussion in Ref. JhJ). A surface insulator-to- metal transition has been 
predicted for the ferromagnetic insulator EuO JTl[: For temperatures below 
Tq = 69 K and decreasing, the majority bd conduction band shifts towards 
the occupied 4/-f bands thereby reducing the insulating gap. This so-called 
red shift is transmitted to an unoccupied Tamm state which is predicted to 
split off at the (100) surface from the lower edge of the conduction band. The 
energy difference between the surface state and the majority bulk band edge 
together with the T-dependent red shift is just sufficient to bridge the gap. 
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Fig. 1. Bulk and the surface density of states (schematic) in case of a metallic 
surface and an insulating bulk. Ep: Fermi energy. 



Therefore, the surface state should become populated for low T which would 
imply a transition to a (half-) metallic state at the surface. 

In any case the concept of a crystal-induced surface state is based on a 
model of effectively independent electrons. It is the detailed form of the one- 
clcctron potential that determines the energy position (lateral dispersion) of 
a surface state in the band structure. Hence, it is obvious that a surface state 
is not normally pinned to E-p (chemical potential) ; one may conclude that for 
a surface of a band insulator there is no a priori reason for a surface state to 
be partially filled. 



4 Mott Transition 

The question is whether or not a situation is conceivable where the appear- 
ance of the surface state is necessarily connected with the metal-insulator 
transition, i.e. where the surface state necessarily appears at Ep. It is now 
clear that to this end one has to look for an electron-correlation effect. 

Consider electrons in a narrow non-degenerate band interacting with each 
other via an on-site Coulomb repulsion U as described by the Hubbard model 
[ fl2| . Any symmetry- broken phases will be excluded from the discussion. At 
half-filling n = 1 and for strong U, the system is then a paramagnetic insula- 
tor. The k-integrated one-electron excitation spectrum (DOS) of this so-called 
Mott-Hubbard insulator |3 13 has the same form as shown in the lower part 



of Fig. [I] where now the two peaks have to be interpreted as the lower and 
the upper Hubbard band separated by an energy of the order of U. 

On the other hand, for U = and for the weakly interacting case, the 
system is a normal Fermi liquid. At an intermediate interaction strength U c 
of the order of the width W of the non-interacting DOS one thus expects 
a metal-insulator transition. This Mott transition is a prime example for 
a quantum phase transition at T — which results from the competition 
between the electrons' kinetic energy ~ W which tends to delocalize the 
electrons and their potential energy ~ U which tends to localize them. 
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The Mott transition is a true many-body effect that cannot be explained 
by a simple perturbational approach. Even within the framework of simplified 
model Hamiltonians, as the Hubbard model, an ultimate theory of the Mott 
transition is still missing ||]. A decisive step forward, however, has been 
made in the last decade with the development of dynamical mean-field theory 
(DMFT) and its application to the Mott transition (see Sec. p§, for a review 
see Ref. p4| , recent results can be found in Ref. Jl5|). 

Within the DMFT one finds that the transition at U = U c is characterized 
by a diverging effective mass m* — > oo or a vanishing quasi-particle weight 
z oc (to*) -1 — > 0, respectively. For U < U c but close to the critical point, 
the DOS has a three-peak structure, consisting of the two well-developed 
Hubbard bands as well as a narrow quasi-particle resonance at the Fermi 
energy with weight z - the DOS has the same form as the "surface DOS" 
shown by the upper part of Fig. |l|. 

It is conceivable that Fig. [[] describes a situation where the bulk of the 
system is a Mott-Hubbard insulator while the surface is in a metallic Fermi- 
liquid state. The quasi-particle resonance would then be a surface state (one- 
electron surface excitation) with a layer-dependent weight z a decreasing ex- 
ponentially with increasing distance from the surface. This surface state nec- 
essarily appears at the Fermi energy as it corresponds to low-energy ex- 
citations well known from quantum impurity systems (Kondo effect). The 
question is for which circumstances this coexistence of the Mott-Hubbard 
insulator and the metallic Fermi liquid can be realized. 

5 Mean-Field Approach 

As it is by no means obvious how to construct a (continuum) Landau theory 
for this problem, the method of choice is to formulate and evaluate a mean- 
field theory for an appropriate discrete lattice model. While for a magnetic 
phase transition one can resort to effective spin models such as the Ising 
model without any detailed knowledge of the electronic structure, the Hub- 
bard model as the minimum model to describe the Mott transition includes 
spin as well as charge degrees of freedom and is thus much more complicated. 
Likewise it is much more complicated to find a proper mean-field theory. For 
example, Hartree-Fock theory, weak- and strong-coupling perturbational ap- 
proaches or decoupling approximations for the Hubbard model are clearly 
inferior compared with the Weiss mean-field theory for the Ising model. The 
latter is non-perturbative, thermodynamically consistent and free of unphysi- 
cal singularities in the entire parameter space. Since the Weiss theory becomes 
exact in the non-trivial limit of infinite spatial dimensions D |l6j , this may 
serve as a simple and precise characterization of what is a proper mean-field 
theory in general. One may therefore hope that the same limit will lead to a 
powerful mean-field approach in the case of the Hubbard model, too. 
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Fig. 2. Left: DMFT self-consistency cycle for the Hubbard model on a semi-infinite 
lattice with layer index a = 1, d (d — * oo). See text for discussion. Right: Layer- 
dependent quasi-particle weight z a for a simple-cubic film of d — 17 layers with 
(100) surfaces. DMFT results for uniform nearest-neighbor hopping t = 1, half- 
filling n = 1, T = and U as indicated. 



That the D — oo limit for a lattice fermion model is well-defined and 
non-trivial in fact, has been proven in the pioneering work of Metzner and 
Vollhardt |r7| . Crucial is a proper scaling of the hopping t oc 1/\D to keep 
the dynamic balance between kinetic and potential energy. To convert the 
abstract definition of a proper mean-field theory into a useful concept for 
practical calculations, it has been important to realize that the D = oo 
Hubbard model can be mapped onto the single-impurity Anderson model 
(SIAM) as now one can profit from various methods available for impurity 
problems. This observation has been made by Georges, Kotliar and Jarrell 
[ ^8[|l9|] . The mapping is a self-consistent one which means that the parameters 
of SIAM depend on the one-particle Green function of the Hubbard model. 

To study surface effects one has to consider a variant of the original Hub- 
bard model. Using standard notations the Hamiltonian reads: 

H= ii||d!,j||/3 c l|| Q CT C j||/3o- + "2" n i\\a<y n i\\a-<y ■ (1) 

iliac 

Here iii labels the sites within a layer a parallel to the surface, a = 1 cor- 
responds to the top surface layer. The model ([[]) differs from the original 
Hubbard model by the mere existence of the surface: The sole effect of the 
surface is to terminate the bulk. For any numerical calculation one has to 
assume a finite number of layers, a = 1, ...,d, i.e. a film geometry, and to 
check the convergence of the results for large d — > oo . 

The generalization of DMFT for surface geometries has been developed 
by Potthoff and Nolting pG]. Again, the limit D — oo may serve as a guide to 
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construct a powerful mean-field theory. Assuming uniform hopping param- 
eters, ti, ]a ,j»p = t between nearest neighbors (i\\,a) and (jy,/9), and using 
the same scaling t oc l/vZ>, the model itself as well as surface effects re- 
main non-trivial for D — > oo. As the different layers parallel to the surface 
must be treated as being inequivalent, the mapping procedure, however, be- 
comes more complicated (see Fig. ^, left). The original many-body problem 
for a semi-infinite lattice with layers a — 1,2, ...,d (with d — ► oo) is self- 
consistently mapped onto a set of impurity problems labeled by the same 
index a = 1,2, ...,d. Each S1AM can be treated independently to calculate 
the impurity self-energy S a (u>). There is, however, an indirect coupling which 
is mediated by the self-consistency cycle: Via the Dyson equation of the lat- 
tice model, the on-site Green function G a (uj) for a given layer a depends on 
all layer-dependent self-energies. The free Green function Q a (oj) of the orth 
SIAM which determines its one-particle parameters is then obtained from the 
DMFT self-consistency condition: Qaifjj)^ 1 = G q ((jj) _1 + S a (uj). 

Fig. H (right) shows the quasi-particle weight z a = (1 — Z^,(0)) -1 as 
obtained from the DMFT using a standard ( "exact diagonalization" ) method 
i jjl) to treat the different impurity problems. The critical interaction U c (d) 
of the d = 17 layer sc(100) film lies close to the bulk critical interaction 
C^cbuik ~ 16 (W = 12 is the width of the free bulk DOS). In the metallic 
phase for U < U c (d) there is a quasi-particle resonance in the interacting local 
density of states for each layer with a finite weight z a . As can be seen in the 
figure, the weight has a pronounced layer dependence. While for small U the 
profile has an oscillating character, it becomes monotonous for interaction 
strengths close to the transition. This is a typical result which is observed for 
films with different surface geometries and indicates a universal behavior of 
the critical profile for U — > U c (d). For both, (7 = 6 and U = 12, the surface 
quasi-particle weight z\ is considerably lower than the bulk quasi-particle 
weight z at the film center. This result is plausible: Due to the reduced 
coordination number at the surface q\ < q, the variance At — q\t 2 of the 
surface-layer DOS is reduced which implies the "effective" interaction U /s/Ai 
to be stronger at the surface compared with the bulk. In this respect the 
surface is "closer" to the insulating phase. Yet, for U — > U c (d) all z a vanish 
simultaneously and there is no surface transition. 

6 Critical Regime 

For systematic investigations of films with different (large) thicknesses, with 
different surface geometries and for different model parameters, a numeri- 
cally exact evaluation of the DMFT requires an effort which is out of scale. 
Fortunately, a simplified treatment of the mean-field equations is possible at 
T = for parameters close to the critical point as has been pointed out by 
Bulla and Potthoff ||. This "linearized DMFT" (L-DMFT) is based on two 
plausible assumptions for the critical regime U — > U c : (i) The effect of the two 
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Fig. 3. Phase diagram for the Mott transition at the (100) and (110) surface of the 
sc lattice, tu (> t = 1): modified hopping within the surface layer a — 1. I7 c ,buik: 
bulk critical interaction. 



Hubbard bands on the quasi-particle resonance can be disregarded and the 
resonance basically reproduces itself in the DMFT self-consistency cycle, (ii) 
The resonance has no internal structure and can be described by a one-pole 
approximation. Although these assumptions are approximate, the L-DMFT 
has successfully passed a number of tests which have been performed by com- 
paring with the full theory and which show that the L-DMFT is well qualified 
to give quantitative estimates for critical interactions and critical profiles as 
well as the correct topology of phase diagrams J22|,^0|,|3) . 

Within the L-DMFT the mean-field equations reduce to algebraic equa- 
tions for z a which involve the electronic model parameters and the system 
geometry. For example, in the case of a surface geometry with qu nearest 
neighbors within a layer, q± nearest neighbors within each of the adjacent 
layers and with uniform hopping t and interaction U except for the hopping 
t\\ 7^ t within the surface layer a = 1, the mean-field equations read: 

36 36i 2 
zi = jp{l\\ hi z i + Q± t 2 z 2) , z a = -jp-[q\\ z a + q± z a+ i + q± z a -i) 

(2) 

for a — 2,3, ...,oo. Equations of this type can be solved analytically or by 
simple numerical means. 

Fig. ^ shows the phase diagram for the Mott transition at two low-index 
surfaces of a semi-infinite D = 3 simple cubic lattice as obtained within 
the L-DMFT. For the unperturbed surface (in = t) there is a single metal- 
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Fig. 4. Surface states (lines) in the low-energy part of the one-electron excitation 
spectrum and bulk continuum (grey region) in the D = 2 Brillouin zone for the 
unperturbed (100) surface (tu = i) and different U — + U c , bulk- Note that the energy 
scale is normalized with respect to the bulk quasi-particle weight z with z —* for 
U — » t7c, buik- 



insulator transition at the bulk critical interaction ("ordinary transition"). 
Within the L-DMFT U c ,bu\k = 6t yjq\\ + 2q±^ = 6\/6- An enhancement of the 
surface hopping tu > t which exceeds a certain critical value, leads to two 
critical interactions: At U — l^buik there is the transition of the bulk of the 
system irrespective of the state of the surface ("extraordinary transition"). 
The surface undergoes its own phase transition to the insulating state at a 
second critical interaction C/ CiSur f > C/c,burk ("surface transition"). The critical 
perturbation c depends on the surface geometry. Multi-critical behavior 
is found for tu — tn jC and U = t/ c ,surf = £4, bulk ("special transition"). 

For £/ Ci buik < U < [/c, SU rf a metallic surface is coexisting with an insulating 
bulk. As is demonstrated by Fig. [ljthis implies the existence of a surface state. 
Fig. ^ shows that even for the unperturbed surface (tu — t) in the metallic 
phase close to the ordinary transition there is a surface state and in the 
limit U — > t/c.buik even an infinite number of surface states. These are split 
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off from the bulk continuum of coherent low-energy excitations the width of 
which vanishes proportional to z oc L/ C ,buik — U for U — > f/ Cj bulk' The critical 
profile of the quasi-particle weight in this limit, z a cx a, causes a strong 
surface perturbation of the low-energy electronic structure which drives the 
surface states. These surface states should be considered as a new kind of 
correlation-induced surface states as for the corresponding non-interacting 
model (with t\\ — t) the occurrence of surface states is impossible. 

Concluding, one finds a phase diagram with the same topology as pre- 
dicted by the Landau theory of surface phase transitions. This analogy can 
even be made quantitative: Consider Weiss mean-field theory for ferromag- 
netic order in the Ising model on a semi-infinite lattice. For T — * To, the 
mean-field equation can be linearized, and one has m a = (J/2T)(q^ m a + 
q±m a+ i + q±m a ^i) which by comparing with Eq. (||) immediately yields 
the following correspondences: m a z a , J/2 36t 2 and T U 2 where J 
is the coupling constant and m a the layer-dependent magnetization. It is ob- 
vious that this analogy has a number of implications for the Mott transition. 

For example, for a thermodynamic phase transition (analogously for the 
T = Mott transition) there are two critical exponents that merely involve 
the critical temperature (interaction strength), the "shift exponent" X s and 
the "crossover exponent" <fi [§,1). They describe the trend of T c (U c ) for films 
with thickness d i— > oo and the trend of T C!Sur f (J7 CiSur f) for the semi-infinite 
system near the special transition, respectively. Within the Laudau theory 
(linearized DMFT) one finds A s = 2 and 4> = 1/2. 

7 Conclusions 

Among different types of metal-insulator transitions at a surface of a single 
crystal, the correlation-driven Mott transition from a normal metal to a para- 
magnetic insulator is distinguished as it offers a comparatively simple route 
to a metallic surface phase of an insulating bulk. Formally, this is expressed 
by the equivalence between the respective linearized mean-field approaches 
to the Mott quantum-phase transition and to the thermodynamic (magnetic) 
phase transition. It should be stressed that the equivalence implies that all 
results of the Landau theory of surface phase transitions have a unique coun- 
terpart for the surface Mott transition. This includes phase diagrams, critical 
profiles of the quasi-particle weight, critical exponents and other critical be- 
havior. In this way a comprehensive and consistent mean-field picture of the 
characteristics of the surface Mott transition is obtained. 

The author would like to thank R. Bulla and W. Nolting for discus- 
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